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Abstract: In this note we show how Dijkgraaf and Vafa's hypothesis relating the exact su- 
perpotential of an A/" = 1 theory to a matrix model can be used to describe all the massive 
vacua of the M = 1*, or mass deformed Af = 4, theory including the Higgs vacuum. The matrix 
model computation of the superpotential for each massive vacuum independently yields a mod- 
ular function of the associated effective coupling in that vacuum which agrees with previously 
derived results up to a vacuum-independent additive constant. The results in the different 
massive vacua can be related by the action of SL(2, Z) on the J\f = 4 coupling, thus providing 
evidence for modular invariance of the underlying A/" = 4 theory. 



1. Introduction 



According to Dijkgraaf and Vafa's (DV) proposal [1-4] , the effective superpotentials for a large 
class of A/" = 1 supersymmetric field theories at finite N are computed from related large- A'" 
matrix models. Some of the best tests of this hypothesis have been for certain mass deformations 
of A/" = 4 supersymmetric gauge theory which preserve A/" = 1 supersymmetry — the so-called 
J\f = 1* theory. The latter theory is obtained by adding masses for all three adjoint chiral 
superfields of the J\f = 4 theory. In addition, one can allow arbitrary polynomial deformations 
of the superpotential involving one of the three adjoint superfields. In this note, we will refer 
to this entire class of theories as jV = 1* theories. These theories provide some of the best tests 
of the DV proposal, since the resulting effective superpotential depends in a non-trivial and 
characteristic way on the coupling constant of the A/" = 4 theory. Up till now the checks have 
involved calculating the effective superpotential in the confining vacuum for the quadratic Tr <l>^ 
deformation [4] and then for arbitrary polynomial deformations and holomorphic condensates 
in the confining vacuum [5] . In the present note we extend this analysis to cover all the massive 
vacua of the theory, including the Higgs vacuum. 

We find that the result of the Tr$^ deformation is an effective superpotential in each 
massive vacuum of the Af — 1* theory computed using the DV prescription from the matrix 
model of the form 



where p is a divisor of N. Once the vacuum- independent constant N'^E2{t)/V2 is added, the 
results in different vacua are related to each other by the action of the modular group, SL(2, Z) 
on r, the J\f — 4 coupling, and the results agree perfectly with earlier computations using 
different approaches [6,8]. What is remarkable about the derivation from the matrix model is 
that the modular properties emerge as a result of the computation without being assumed at 
the beginning. 

Let us first briefly recall the vacuum structure of the M — 1* theory. In J\f — 1 language, 
the M — 4 SU(A^) theory with couphng constant r = 4m/gYM + ^/Stt has three adjoint chiral 
superfields $~ and We consider a general class of deformations of the M — 4 theory 
specified by a tree level superpotential 
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where V{^) is a general polynomial 
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Note that for the sake of simphcity we have set all masses to unity in the knowledge that it is 
simple to re-introduce them. For V{^) = $^ we obtain the basic A/" = 1* theory, however, it is 
also interesting to consider the space of possible J\f = 1* deformations as in [5]. 

The vacuum structure can be deduced by proceeding with a classical argument augmented 
by some quantum considerations. The classical vacuum structure follows by solving the F- and 
D-flatness conditions modulo gauge transformations. Equivalently, one solves the F-flatness 
conditions modulo complex gauge transformations. Up to complex gauge transformations, the 
solutions of the F-flatness conditions are associated to representations of SU(2) in the following 
way: $ is precisely iJ^ of an SU(2) representation, reducible or irreducible, of dimension N, 
and $^ have the same non-zero elements as J* except that the actual numerical values depend 
on the form of the potential V{^). Vacua can be classified according to whether there are, after 
the Higgs mechanism, unbroken U(l)'s, in which case it is a massless vacuum. On the contrary, 
if the unbroken gauge group is non-abelian (or empty) one has a massive vacuum. In terms 
of SU(2) representations, the massive vacua correspond to N/p copies of the p-dimensional 
representation, where, by construction, p must be a divisor of N. The unbroken gauge group 
is then SU {N/p). The classical eigenvalues of $ are explicitly 

Xf = i(p-2j + l) , j = l,...,p, (1.4) 

each with a degeneracy of N/p. Now quantum reasoning is needed to deduce the multiplicity 
of vacua for each p. One expects at low energies that the theory flows to pure A/" = 1 super- 
symmetric gauge theory with gauge group SU (N/p). Standard arguments then imply there are 
N/p physically inequivalent quantum vacua. All-in-all, there are "^p^j^N/p = J2p\nP massive 
vacua, where p is a divisor of A^. 

It is the goal of this paper to calculate ( p..l|) from the matrix model proposal of Dijkgraaf 
and Vafa [1-4], generalizing the results of [4,5], which held for the confining vacuum {p = 1 
above). 

2. The matrix model and its solution 

As a direct consequence of the proposal of Dijkgraaf and Vafa [1-4], in a given vacuum the exact 
effective superpotential for the above class of deformations of the A/" = 4 theory is computed 
by the planar diagram expansion, i.e. a large- A^ limit, of the matrix model for which there is a 
matrix variable for each chiral multiplet of the field theory and whose action is the tree-level 
superpotential of the field theory. Using the same notation for the matrix variables and their 
associated superfields, the matrix model partition function is 

Z= /"[d$+][ci$-][d$]exp-— Tr(i<l>[<I)+, $-] + $+$- + r($)) (2.1) 
J 9s 
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In the matrix model, unlike in the field theory, one takes <l>^ = ($^)^ and the fiuctuations of 
$ around the ensuing saddle points to be Hermitian. At the saddle-point we will allow the 
eigenvalues of $ to be complex. As noted for example in [10], the key fact that permits the 
solution of these matrix models is that one may integrate out $^ exactly to obtain a one-matrix 
integral to be solved in the large- limit: 

Z= [d^-^—-, -. (2.2) 

Met(Adj$ + ^) ^ ^ 

The resulting one-matrix model actually becomes tractable in the large-A^ limit by going to 
the eigenvalue basis and performing a large- A^ saddle-point approximation to the integral. The 
details of this procedure have been extensively discussed in the literature. We refer the reader 
to [10], and references therein, for more details. 

In order to describe the confining vacuum, for which $ = classically, one takes a one cut 
solution as follows [4, 5] . The eigenvalues Xj of $ interact via a repulsive effective potential 
and form a continuum in the large- A^ limit and condense onto a cut along the real axis. The 
actual extent of the cut and the density of eigenvalues p(A) along the cut is self-consistently 
determined by the saddle-point equation in terms of the parameters of the deformation V{^) 
and the matrix model 't Hooft couphng S = gsN. Proceeding along these lines leads to a 
description of the confining vacuum [4,5]. 

For a general massive vacuum associated in the nomenclature of SU(2) representations to 
N/p copies of the p-dimensional representation, the classical eigenvalues of $ are given in (|1.4| ). 
This suggests in the matrix model we should take a multi-cut solution where at large A^ the 
eigenvalues condense on p cuts defined as 

Qj = = \f + X , -aj < X < I , j = l,...,p . (2.3) 



The saddle-point equation in the large- A^ limit is most conveniently written in terms of the 
resolvent function 

u;{z) = [ -^dA (2.4) 

J Z — A 

where p(A) is the unit normalized spectral density of eigenvalues: 

1 ^ 

i=l 

a function that in the large-A^ limit only has non-trivial support on the cuts. The spectral 
density is normalized in such a way that the filling fractions of eigenvalues along each cut are 
given by 

i\T. r<^j 

p{Xf + x)dx. (2.6) 
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where Nj denotes the number of eigenvalues along the jth cut. 



The resolvent uj{z) is an analytic function on the complex 2;-plane which has cuts precisely 
along each Cj, j = 1, . . . ,p. The discontinuity across each cut gives the spectral density at that 
point: 

uj{X + ie) -uj{X-ie) = -2mp{Xy, X e Gj . (2.7) 



In terms of uj^z), the saddle-point equation is 
V'{X) 



uj{X + ie) + uj{X — it) 
— I {uj{X + i + ie) + uj{X + i — ie) + uj{X — i + ie) + uj{X — i — ie)) 



(2.8) 



for A G UjCj. Notice that various principal values have been taken, giving rise to the ±ie 
prescriptions, in order to render the integrals well defined. Unlike the single cut solution, the 
final terms — those shifted by ±i — also have to be defined as a principal value since most of the 
shifted cuts collide with each other in pairs. 

At this point we need to employ some extra guesswork. We start with the observation 
that in the massive vacuum there should be same number of eigenvalues on each cut. In other 



words the filling fractions (2^) should be the same. This is well motivated by the classical 
solutions for massive vacua of the field theory wherein the adjoint scalars have VEVs that 
split into N/p copies of a p-dimensional representation of SU(2) preserving an SU{N/p) gauge 
symmetry classically. Of course, we would also expect solutions of the saddle-point equations 
with different filling fractions which would correspond to massless vacua. However, we are going 
to specialize to the massive vacua by making an ansatz which automatically ensures that the 
filling fractions are identical. The ansatz we make is that the eigenvalue density is the same 
along each of the cuts (all of which have the same extent aj = a): for each pair j, k 

p[Xf + x) = p{Xf + x) , -a<x<a, (2.9) 

or equivalently 

uj{Xf + X + le) - uj{Xf + x-ie) = uj{Xf + x + ie) - u;(A^' + x - ie) . (2.10) 



Having made this ansatz we now investigate the saddle-point equations (p^.8|). Notice that 
they are linear in the w{X + ie) + w{X — ie) and so we can solve for these quantities at each cut. 
One finds 

u{Xf + x + ie)+ uj{Xf + x-ie) 

= 2Y{Xf + x) + ^ ((p + 1 - j)uj{x + i{p + l)/2) + Mx - i{p + l)/2)) , ^^'^^^ 
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where —a < x < a and j = l,...,p. Notice that the resolvent has no discontinuities at 
x±i{p+ l)/2. In the above, Y{X) is the polynomial which is determined by the potential ^(A) 
via 

= 2r(A) - r (A + - r(A - 1) . (2.12) 

As in the one-cut solution, it useful to define the function 

G{z) = U{z) + iS{u{z + i) - uj{z - f )) (2.13) 
where U{z) is a polynomial in z such that 

iV'{z) = U{z + ^)-U{z-l) . (2.14) 



Given the ansatz ( |2.9| ), the analytic structure of G{z) is much simpler than that of the resolvent 
itself. The cuts of uj{z) are shifted up and down by ±| and because of ( p.lOj ) most cancel 
pairwise. There are only two remaining cuts — the ones on the outside of the stack — given by 
+ X and — |p + X, respectively, for —a < x < a. Furthermore, putting ( p.ll| ) and ( p.lO|) 
together, we end up with a rather simple and suggestive constraint on G{z): 

(^(x + |p ± ie) = G(x - |p =F ie); -a < x < a . (2.15) 



Now we see that having made the ansatz we now have what is closely related to the one-cut 
solution, p = 1, of [4,5]. In terms of the function G{z), the only difference is that the two cuts 
of G{z) are shifted from ±i/2 to ±ip/2. Furthermore, just as in the p = 1 case, the equation 
( |2.15| ) has the effect of gluing the top (bottom) of the upper cut to the bottom (top) of the 
lower cut making what is effectively a torus. This is very important, since on general grounds 
we expect each of the massive vacua to be associated to a particular torus with its own complex 
structure [6,7,11,12]. 

Fortunately, we can now almost lift the analysis of the one-cut solution, with only slight 
modifications to account for the new positions of the cuts, and we refer to [5] for more details 
and conventions. Firstly, we can uniformize the torus in the w-plane by means of the mapping 



z{u) = 



TT 



C{u) U 

UJi 



(2.16) 



The torus is defined by the quotient of the complex u-plane by the lattice generated by the two 
complex numbers 2uJi and 2uj2. The complex structure of the torus f = 002/ uji is a parameter. 
In effect we shall see that we have changed variables from a — and hence from S — to the new 
variable f. Notice in ( p.l6| ), the multiplicative factor of p which ensures that the cuts occur in 
the requisite positions. Tightly encircling around the upper (lower) cut anti-clockwise in the 
2;-plane corresponds to LJ2 + 2xu!i (or —U2 + 2(1 — x)u!i) for < x < 1 on the w-plane. The 
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former, or minus the latter, defines the A-cycle of the torus. In addition, the point at infinity 
in the z plane corresponds to m = 0. 



At this point we find it convenient to choose the simplest A/" = 1* deformation obtained by 
taking a potential = although it should be said that more complicated potentials are 

just as tractable [5]. In this case, Eq. ( |2.13|) fixes the asymptotic form of G to be 

G{z) ™ z^^O(\lz^) . (2.17) 

Given the fact that G must be single valued on the torus — an elliptic function — this fixes its 
form uniquely. Eq. ( p. 17] ) implies that in the u plane G{z{u)) must have a double pole at m = 
with a set coefficient. Given that ( p.l6|) implies 

«=-f^^i + PMf^i + ..., (2.18) 

7t Z TX-^ Z-^ 



means that 



G(.M) = P^#U.)-^". (2.19) 



Once again take note of the factor of which will prove important. 

We are now ready to implement the Dijkgraaf-Vafa prescription. First of all, recall that 
the filling fractions for each of the p original cuts are given as in Eq. (|2.6|) . For our ansatz ( p. 9]) 
each of the filling fractions are equal and following the analysis in [4,5] we can express them as 
an integral of G{z{u)) around the A-cycle of the torus: 

2mS, = -i f G{z{u))'^^du = ^^^^ , S, = gsN, . (2.20) 

J dU iZ (XT 

The factor of arising from the p^ and p factors of G{z) and z(m), respectively, will turn out 
to be crucial. In the above E2 is the 2^^*^ Eisenstein series (see [5] for definitions). 

The other important ingredient in the Dijkgraaf-Vafa proposal is the quantity dJ-'o/dSj 
being the variation in the genus zero free energy of the matrix model upon transporting an 
eigenvalue in from infinity to the j^^ cut to infinity. In the present context where we will take 
all the filling fractions to be the same, we only need an expression for the sum X]j=i dJ-'o/dSj. 
This quantity is then simply related to an integral of the same form as (|2.2CI|) , but now over the 
B cycle of the torus —ui + (2x — l)uj2, < x < 1. The integral is easily evaluated following [5] 

t i OiMu)/-^^U . .3(^^ _ J_,,,)) . ,,,, 

j = l J 

Again the factors of p are crucial: G{z) and z{u) give p^ and p respectively. 
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The effective superpotential is obtained by extremizing the following expression with respect 



We^ = ^2{N,^,-2^^rS^ . (2.22) 



to Sf 

p 

In this expression Nj are the number of eigenvalues of $ which are classically equal to A^^ It is 
important to realize that these are not directly related to the A'^ in the matrix model (although 
we use the same notation). The latter are related to the filling fractions Sj/ S via Nj = Sj/gg- 
In order to describe the massive vacua we now take the physical Nj = N/p. Since Nj is an 
integer, p must be a divisor of A^. In this case, after having made the ansatz (|2.9|) , which 
solves the saddle-point equations, each of Sj depends upon only one parameter, namely f, the 
complex structure of the auxiliary torus in the w-plane. Hence, in implementing the Dijkgraaf- 
Vafa prescription for the massive vacua, we only have to extremize the superpotential with 
respect to f. Substituting our expressions ( p. 201 ) and ( |2.21| ), we have 



Wes = Np" (^)^^^ - 2mTp^S,=i . (2.23) 

where we have written the result in terms terms of the (p = 1) one-cut solution quantities in 
order to highlight the p dependence: 

^ 12 df ' V dS Jp=i 12 df 12 ^ ^ ^ ' 



It is easy to see that 

and substituting back we have 



dWes ^ ^ p^T 

^ = ° ^ ^= AT 



:2.25) 



Weff = -^E2{p't/N) . (2.26) 



where as stated above p has to be a divisor of A^. Upon restoring the mass scales in the problem 
this result agrees precisely, up to a vacuum- independent additive constant, with the expressions 
of [6,7]. Note that SL(2, Z) modular transformations of r relate the superpotentials in different 
massive vacua with the understanding that Wes has modular weight 2 and that one must 
account for the anomalous modular transformation of E2 by adding the vacuum independent 
constant N'^E2{t)/12} This is remarkable in that ^-duality of the underlying Af = 4 theory 
was not assumed to begin with, but rather emerges from the solution of the matrix model. 

We want to emphasize that our method also covers the Higgs vacuum (and other vacua 
where p scales with A^) which at first sight is puzzling because in the field theory all the eigen- 
values are — classically at least — non-degenerate. How can they form continua on A^ separate 

^The freedom to add such a constant is an example of the more general effect of operator mixing (in this 
case with the identity) [8] . 
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cuts in the matrix model when there are only eigenvalues? The resolution of this puzzle is 
the realization that the Nj in the field theory are not equated to the Nj = Sj/qs in the matrix 
model. The latter can be large while the former can be set to one in (|2.22|) . 



Notice that our result ( |2.26| ) only describes one of the N/p massive vacua associated in 
the original analysis to the p-dimensional representation of SU(2). The other N/p — 1 such 
vacua are obtained by repeatedly performing the modular transformation r ^ r + 1 a sufficient 
number of times to give Ql.lj) . 

Our multi-cut technology can also be employed in the Af = 1* deformation of the Leigh- 
Strassler theory for which the confining vacuum result was quoted in [5] . The result in the p*^ 
massive vacuum is 

^ N cos(/3/2) e,i(3p/2\p'T/N) N 

p4sin3(/3/2) e[{(3p/2\p^T/N) 4sm\(3/2) ' ^' ' 

Acknowledgments S.P.K. would like to acknowledge support from a PPARC Advanced 
Fellowship. 



References 

[1] R. Dijkgraaf and C. Vafa, arXiv:hep-th/0206255. 

[2] R. Dijkgraaf and C. Vafa, arXiv:hep-th/0207106. 

[3] R. Dijkgraaf, Talk at "Strings 2002" , Cambridge, 

http : / / www . damtp . cam . ac . uk/ str ings02/ avt/di j kgraaf / 

[4] R. Dijkgraaf and C. Vafa, [hep-th/0208048] . 

[5] N. Dorey, T.J. Hollowood, S. Prem Kumar and A. Sinkovics, [hep-th/0209089] . 

[6] N. Dorey, JHEP 9907, 021 (1999) [hep-th/9906011] . 

[7] N. Dorey and S. P. Kumar, JHEP 0002, 006 (2000) [hep-th/0001103] . 

[8] O. Aharony, N. Dorey and S. P. Kumar, JHEP 0006, 026 (2000) [hep-th/0006008] . 

[9] N. Dorey V. V. Khoze and M. P. Mattis, Phys. Lett. B 396, 141 (1997) [hep-th/9612231] . 

[10] V. A. Kazakov, I. K. Rostov and N. A. Nekrasov, Nucl. Phys. B 557, 413 (1999) 
[hep-th/9810035] . 

[11] N. Dorey T. J. Hollowood and S. Prem Kumar, Nucl. Phys. B 624, 95 (2002) [hep-th/0108221] . 
[12] N. Dorey and A. Sinkovics, JHEP 0207, 032 (2002) [hep-th/0205151] . 



8 



